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We derive exact expressions for the local entanglement entropy £ in the ground state of the one- 
dimensional Hubbard model at a quantum phase transition driven by a change in magnetic field 
h or chemical potential (i. The leading divergences of d£/dh and d£/dp are shown to be directly 
related to those of the zero-temperature spin and charge susceptibilities. Logarithmic corrections 
to scaling signal a change in the number of local states accessible to the system as it undergoes the 
transition. 
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Entanglement is a generic feature of quantum systems, 
implying the possible existence of non-local correlations. 
Such correlations — which lead to highly counter-intuitive 
phenomena — were long seen as an artifact of quantum 
mechanics |jj . With the advent of quantum information 
theory it is now understood that entanglement — and 
the correlations associated with it — is not only intrinsic 
to the fabric of reality 0, but can also be used as a 
physical resource, essential for performing such tasks as 
teleportation or quantum computing . 

A new line of research 0, points to a connection 
between the entanglement of a many-particle system — 
as quantified by a properly chosen measure — and the 
appearance of a (zero-temperature) quantum phase tran- 
sition (QPT) 0. Barring accidental occurrences of non- 
analyticity, a discontinuity [singularity] in the [derivative 
of the] ground state concurrence of an A^-qubit system 
appears to be associated with a first [second] order QPT 
[jj (with concurrence measuring the entanglement be- 
tween two neighboring qubits These and related 
results are important as they hold promise of novel per- 
spectives on condensed matter, drawing on insights from 
quantum information theory. By analyzing entanglement 
properties one expects to gain insight into how the as- 
sociated non-local (purely quantum) correlations influ- 
ence the critical behavior of a quantum phase transition. 
Building an understanding of this connection should en- 
able breakthroughs in the design of future experimental 
probes of collective quantum phenomena. Also, archi- 
tectures for quantum information processing that take 
advantage of the entanglement in the vicinity of a quan- 
tum phase transition (quantum adiabatic computing) 
should benefit from a detailed understanding of entan- 
glement scaling properties. 

Most results to date on the entanglement - QPT con- 
nection have been obtained from numerical studies of 
finite lattice s pin systems, supplemented by some ana- 
lytical results [l(j]. Much less is known about entangle- 
ment scaling properties of itinerant electron systems. In 
this Letter we make a dent on this important problem 
by studying the one-dimensional Hubbard model close 
to a quantum phase transition. Recent work on this 
and related models show that features of the ground 



state phase diagram can be reproduced by studying cer- 
tain characteristics of the local entanglement entropy 
[H EIIl El 03 Here we exploit the Bethe Ansatz 
solvability of the Hubbard model to derive exact expres- 
sions for the critical scaling of the local entanglement 
entropy £(ipo) of its ground state \ipo) as function of mag- 
netic field h and chemical potential /i. We find that the 
leading scaling behavior of <9£/<9/i for repulsive interac- 
tion coincides with that of the charge susceptibility xc- 
A similar result holds for d£/dh, but with logarithmic 
corrections that signal a change of dimension of the ac- 
cessible local state space at the transition. The fact that 
an entanglement measure of a critical many-particle sys- 
tem can be quantitatively linked to a physical observable 
is a striking result, and goes beyond standard construc- 
tions of entanglement witnesses 0] that merely detect 
the presence of entanglement. To what extent our results 
can be generalized to other quantum systems is yet to be 
answered. 

To set the stage, let us recall that the very notion of en- 
tanglement of a composite quantum system relies on the 
tensor product structure of its Hilbert space. When the 
system is made up of itinerant electrons, however, the 
physical subspace is restricted to an anti-symmctrized 
one which lacks a natural product structure. One may 
circumvent the problem by passing to an occupation 
number representation spanned by the 4 L basis states 
I h)i<8 I ri)2 ® — <8 | ti)l, where, in obvious notation, 
\ n )j H lT)ii IDji or IT4)j is a local state at site j, 
with L the number of sites on the lattice ^3- This is 
a convenient basis in which the tensor product struc- 
ture is manifestly recovered, with the local states describ- 
ing electronic modes easily accessible to an observer. By 
splitting the system into two parts A and B one can then 
proceed as usual and define the entanglement entropy £ 
of a pure state | ip) as £ — — Tr(p^log 2 pA) 0- The re- 
duced density matrix pa = Tr^(p) is obtained from the 
full density matrix p = \ ip}{ip | by tracing out the local 
states belonging to B. In what follows we focus on the 
entanglement entropy of a single site, obtained by taking 
A to be a single (arbitrarily chosen) site, with B the rest 
of the system. 

We begin by studying the Hubbard model with an ap- 



2 



plied magnetic field H: 
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H = - t E C ]aC j+ S a + Uj2 n rt n n-»BHj2S?. (1) 
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Here cj a and c Q are creation and annihilation operators 

for electrons at site j with spin a =1 , J,, and nj a = c ] a c j a 
and 5| = (n-jf — nji)/2 are the corresponding num- 
ber and spin operators, respectively. We assume peri- 
odic boundary conditions. In the following we shall work 
with dimensionless quantities u = U /At and h = fisH/t, 
putting t = 1 (of dimension energy). Using the fact that 
the Hamiltonian in Q is translational invariant and con- 
serves particle number as well as the z-component of the 
total spin, it is easy to verify that the reduced density 
matrix pa for a single site A is diagonal in the chosen 
basis. It follows that the corresponding single-site en- 
tanglement of the ground state | tpa) can be written as 



£ = -Wq l0g 2 U>0 — W| log 2 Wf — Wl l0g 2 U>1 — W 2 l0g 2 W 2 

(2) 

with 

w 2 = (nj^rij^o, w a = (n ja ) - w 2 , a =f, J. 

Wq = 1 — lOf — Wl — W 2 ■ (3) 

The problem is thus reduced to calculating the expecta- 
tion values for double and single (spin-up and spin-down) 
occupancies in the ground state. 

Let us first look at the case of attractive interaction, 
u < 0, with n = 1 (half-filling). In the limit \u\ ^> 1 
we can use the Hellman-Feynman theorem, (dH/du)o = 
OEq/Ou, together with the known Bethe Ansatz result 
for the ground state energy (lij. Eq/AL = u{l/2 — m) — 
(l/27r)sin(27rm) +0(l/u) to obtain 



w 2 



li>E " l- m + 0(l/u% (4) 
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with m = (1/2L) 2^=1 ( n iT ~ n jl) ^ nc magnetization per 
site. Neglecting the 0(l/u 2 ) corrections it follows imme- 
diately from J3J and Q that 

W| = 2m, luj, = 0, wo = g — m ' ft > 0. (5) 

Combining ©,0, and JHJ), we obtain for the single-site 
entanglement: 

£ = -2mlog 2 (2m) - (l-2m)log 2 ( m), h>0. (6) 

The dependence of the magnetization on the applied field 
can also be derived from the ground state energy, and one 
finds 

< h< h cl 

i(h) = { ^arccos(-(u+ |)) h cl < h < h c2 (7) 
| h c2 < h 



with lower [upper] critical field h c \ = A(\u\ — 1) [h c2 = 
A(\u\ + 1)]. The single-site entanglement as a function 
of magnetic field, £ — £(h), can now be read off from 
© and 0. The result for the |u| > 1 limit is plotted 
in Fig. for large values of h. Note that in this limit 
there are two local states, | 0) and |t|), available to the 
system when h < h c i, implying that £(h) = 1. In con- 
trast, the fully magnetized state for h > h c2 is a direct 
product of local spin-up states, and hence £ (ft) = 0. For 
comparison we have plotted the single-site entanglement 
for free electrons also in Fig. (for both positive and 
negative values of the magnetic field) . This result is eas- 
ily obtained from Ref. by noting that w 2 — 1/4 — m 2 
when u = 0, with m = (1/ir) arcsin(ft/4) in the interval 
-A<h<A. 




Magnetic field 

FIG. 1: Entanglement entropy £ of a single site versus mag- 
netic field h for the attractive Hubbard model with \u\ 2> 1 
(solid curve). For comparison, the single-site entanglement 
for the free case (u — 0) is shown by the dotted curve (on a 
different scale). 



The phase transitions at h c \ and h c2 are second or- 
der, with diverging spin susceptibilities xsi — (327r 2 |/i — 
ftci|) -1 ^ 2 j i = 1,2 [18j. As mentioned in the introduction, 
in the case of an ./V-qubit system there is strong evidence 
that the derivative of the ground state concurrence with 
respect to a critical parameter diverges at a second-order 
QPT 7J. For the present problem (where the local de- 
grees of freedom have four components, not two as for 
a qubit) the plot in Fig. suggests a divergence of 
d£/dh as h — ► h c \ + and h — * h c2 ^. To analytically check 
whether the derivative of the single-site entanglement is 
a true marker of quantum criticality for this problem, we 
write u + ft/4 = (h — h e i)/A + (— l) 1 , i = 1, 2 and expand 
d£ / dh in h — h c \ and h c2 — h, respectively. We obtain 



XSi 



- = (-D 
dh 1 ' ln(2) 



(In | ft - h oi \ + const.), i = l,2 (8) 



for h — * h c i + and h — > ft C 2— j respectively. This confirms 
that d£/dh diverges at the magnetic phase transitions. 
Moreover, it shows that the divergence of d£/dh is given 
by the spin susceptibility — up to a logarithmic correction. 
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By writing the leading term of d£/dh on the "mixed" 
form (—dm / dh){2\a(w^) — ln(u> ) — ln(w 2 ) ) / In 2 [cf . Eqs . 
and 10], and combining this expression with J7J, one 
sees that the logarithmic divergence in JSJ) comes from 
a change of the number of local states accessible to the 
system as it undergoes the transition: As h — ► h c \+ the 
local spin-up states get suppressed (u>j — * 0), while for 
h — ► h c 2- both empty and doubly occupied local states 
get suppressed (wq,W2 — + 0). 

Turning to the half-filled case with repulsive interac- 
tion, u > 0, a QPT now occurs only at the value of 
the field for which the magnetization saturates: h C 2 — 
4(\/l + u 2 — u) |19|. As shown by Takahashi, the ground 
state energy for any finite value of u > in the criti- 
cal region h — ► h C 2- can be expanded in terms of the 
expectation value for single spin-down occupancy |2(i| : 

E, 



7T 2 1 



■.(n 3i )l + 0((n 3i )t). (9) 



24 *JTTu- 

With the same procedure as used for the attractive case 
above, Eq. ©, together with © and J3J), yield: 



d£ C 



dh 21n(2 



■Xs (In \h-h c2 \+ const.), h^h c2 -. (10) 



Here C = 2 - u/Vl + u 2 , and 2ttxs = (4 + 4m 2 ) 1 / 4 !/! - 
/i^l" 1 / 2 - The logarithmic correction in now signals 
the suppression of all but the spin-up states as one ap- 
proaches the saturation point h C 2 from below. 

We next study the effect of a varying chemical potential 
on the single-site entanglement of an open system. We 
make the simplifying assumption that the environment 
acts solely as a particle reservoir j^H, and add the term 



■^E,-=i(%- 
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to the Hamiltonian in QJ , with 



\x a dimensionless chemical potential (multiplied by the 
hopping amplitude t — 1). To simplify further we turn 
off the magnetic field in (JIJ, putting h = 0. Focusing on 
the case of repulsive interaction, u > 0, with n < 1, 
the system exhibits two quantum critical points ~|22|: 
/i c i = — 2 and /i C 2 = 2 — 4 J Q Ji(u))(u[l+eKp(2uu)]) du>, 
with Ji(w) a first-order Bessel function. Both transitions 
are second-order with diverging charge susceptibilities 



XCr = C(u)\fl 



/'<-: 



-1/2 



1, 2 in the limits /i — > /Lt c i+ 



(empty lattice transition) and /x — > (J- C 2- (Mott transi- 
tion), respectively (with c(u) a positive u-dependent con- 
stant). To obtain the single-site entanglement £ we first 
notice that 7i M conserves spin and particle number for 
fixed /i, and that hence the expression for £ in J2J) re- 
mains valid. Recalling from the Lieb-Mattis theorem [2^] 
that the ground state has zero spin (for any n with nL 
an even integer) we can write the parameters appearing 
in J5J) as 



w 



1 — n + W2, w\ = wi 



n 
2 



w 2 - 



(11) 



The value of W2 can again be extracted from the ground 
state energy via the relation W2 = (dEo/du)/4L, where 
the Bethe Ansatz solution for E can now be expressed 
via a 1/u expansion |2^j : 



1=1 



(12) 



The values of «/(n) are tabulated to fifth order in Ref. 
[24} . The ground state energy in (|12(l also determines the 
chemical potential as function of filling: (J,(n) — dEo/dn. 
By inverting fj,(n) and inserting the resulting values for 
the w-parameters from Q into J2J) we can plot £ vs. /i for 
any value of u > 1 . Some representative plots are shown 
in Fig. J5J), together with the single-site entanglement 
for free electrons (u = 0). 




Chemical potential 



FIG. 2: Entanglement entropy £ of a single site versus chemi- 
cal potential n for the repulsive Hubbard model. The plateaus 
correspond to half-filling (n = 1), cut off at [i = 2. The dotted 
curve is that for free electrons (u — 0), plotted in the region 
< n < 2. 



In order to analytically explore the quantum critical 
regions /1 — * /i c i+ and /i — > \x C 2- we first consider the 
u — > 00 limit where = 0. In this limit l|12|) implies that 
n(fi) = (1/ 7r) arccos(— fi/2). Combining this expression 
with Eqs. @ and we obtain 



dfi [ ' 21n(2) 



(In — /z c j| + const.), i — 1,2 (13) 



for /1 — > /i c i+ and // — > /i c 2-, respectively. Note that 
the derivative of the single-site entanglement is again 
given by a susceptibility, corrected by a logarithmic fac- 
tor that reflects the change in the number of available lo- 
cal states as the system undergoes the transition: When 
M ~~ * \p — * A*c2-] the singly occupied [empty] local 

states get suppressed. (Cf. the argument after Eq. JSJ.) 
The exact analogy with the magnetic scaling in (JSJ can 
be understood by carrying out a particlc-holc transfor- 



mation for spin-up electrons: 



(— l) J c T (leaving 



the spin-down electrons untouched). This transforma- 
tion maps the zero-field attractive system with a variable 
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chemical potential onto a half-filled repulsive system with 
a variable magnetic field, ft follows that the single-site 
entanglement at /x cl (/i C 2) has the same behavior as at 
h c2 (h cl ). 

Turning to the case of large but finite u, we focus on the 
Mott transition /i — ► fi C 2-- A straightforward analysis, 
again using the Bethe Ansatz result in (|12|l . yields for the 
leading behavior of the single-site entanglement: 



in the original Fig. 2 of this paper (published in Phys. 
Rev. Lett. 95, 196406 (2005)) |U. H.J. acknowledges 
support from the Swedish Research Council. 



dS 



7T = -C{u)x 



C2, 



(14) 



with C(u) a positive w-dependent constant. Note that 
there is no logarithmic correction in l|14|) : When u is 
finite the metallic (// < and insulating (/i > /i C 2) 
ground states for < n < 1 are both superpositions of 
all four types of local states | 0)j, |t)j, and |T4)j- It 
follows that none of the ui-parameters in (j2J tend to zero, 
and the logarithmic terms add up to a constant |25j . 

The results in lJT5j l and (|T4"|) derived for < n < 1 can 
be extended to the region 1 < n < 2 via the particle- 
hole transformation c\j Q (— l)^Cj- a , a =|, J,. One finds 
that the change of £ at the quantum critical point /i C 3 = 
Au — fi C 2 \p c i = 4it — /i c i] where the system goes to more 
than half-filling, n > 1 [complete filling, n — 2] exhibits 
the same scaling as at the transition to half-filling [empty 
lattice] studied above. 

To summarize, we have found that the derivatives 
d£ I dh and d£ /dfi of the single-site entanglement £ are 
faithful markers for QPTs in the Hubbard model driven 
by a change in magnetic field h or chemical potential 
/i. Via an analysis based on the Bethe Ansatz solution 
of the model, we have derived exact expressions for 
d£/dh and d£/d/j, at criticality, revealing that these 
quantities scale with the corresponding diverging spin 
and charge susceptibilities, respectively. Logarithmic 
corrections signal a change in the number of available 
local states at the QPT. That a critical entanglement 
entropy is directly connected to a susceptibility is an 
intriguing property. A zero-temperature susceptibility is 
an observable determined only by the dependence of the 
ground state energy on the critical parameter (magnetic 
field or chemical potential). The entanglement entropy, 
on the other hand, carries information about the very 
structure of the ground state. From our analysis it is 
transparent how the linkage formally comes about for 
the Hubbard model: The reduced density matrix can be 
parameterized by expectation values derivable from the 
ground state energy via the Hellman-Feynman theorem. 
To find a physical interpretation of the connection, and 
to explore whether it can be extended to other critical 
quantum many-particle systems, is an interesting and 
challenging problem. 

We thank R. Rehammar and O. Syljuasen for valuable 
discussions. We are also grateful to K. Capelle and V. 
V. Franca for drawing our attention to a numerical error 
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